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A number of practical problems, including questions about reliability of 
Geiger counters and short-circuits in electric cables, reduce to the mathe- 
matical problem of coincidences in Poisson patterns. This paper presents 
the probability of no coincidences as well as asymptotic formidas and simple 
bounds for thai probability under a variety of circumstances. The probability 
of exactly N coincidences is also found in some cases. 

INTRODUCTION 

A number of prai^tical problems are questions about what we call 
"coincidences" in Poisson patterns. In rf-dimensional space, a Poisson 
pattern of density X is a random array of points such that each infinitesi- 
mal volume element, dV, has probability \dV of containing a point, 
and such that the numbers of points in disjohit regions are independent 
random variables. Then a volume, V, has probability 

of containing exactly k points. A coincidence, in our usage of the word, 
is defined as follows: We imagine a certain fixed distance 5 to be given 
in advance; two points are then said to be coincident if they he within 
distance 5 of one another. 

Examples 

The best-kno\TO case of a coincidence problem concerns Geiger coun- 
ters. In the simplest mathematical model, there is a short dead-time 5 
after each count dui-ing which other particles can pass through the 
counter \\"ithout registering a count. In our present terminology, a coimt 
is missed whenever two particles traverse the counter with coincident 
times of arrival. The same problem is encountered with telephone call 
registers. 

1005 
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Another example arises in the manufacture of electric cable. Each 
wire in a cable is covered with an insulation which contains occasional 
flaws. When the cable is assembled it will fail a short circuit test if it 
contains a pair of wires such that a flaw on one wire Ues within some 
distance 5 of a flaw on the other wire. In a similar way, coincident flaws 
in the insulation of the wire from which a coil is wound can lead to fail- 
ure of the coil. 

There are also some problems in the development of certain military 
systems which lead to the consideration of coincidences in Poisson pat- 
terns. 

Outline of Work 

Our primary aim is to study the probability of no coincidences under 
various circumstances. In Part I, we examine coincidences of two differ- 
ent Poisson patterns, of densities X and n respectively, on a line of length 
L. Here we do not count two points of the same pattern within a distance 
5 as giving a coincidence. A set of integral equations yields the probabihty 
of no coincidences as well as an asymptotic formula and upper and lower 

bounds. 

In Part II, we study the probability, Fo{L), of no comcidences for a 
single one-dimensional Poisson pattern of density X. These results may 
also be interpreted as the distribution fmiction for the minimum distance 
between pairs of points of a Poisson pattern. Sample formulas are the 
asymptotic formula (for large L) 

FoiL) 



(X - a)[l -H 5(X - a)] 
and the bounds (valid for all L) 

(l - ^) e-^-e-^-^)^ ^ F,iL) ^ e^^^e-'"-''', 

where s = — a is the largest real root of 

The problem of n Poisson patterns, all of the same density X, is ex- 
amined in Part III. Coincidences are now counted between points of any 
two distinct patterns. 

The one-dimensional problems of Parts I-III succumb readily to ana- 
lytic techniques. We can find exact expressions for the probabilities of 
no coincidences in Parts I-III. Two entkely different methods of deriving 
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exact results are available and are illustrated in Parts II and III. Un- 
fortunately, the exact formulas, although they arc finite sums, contain a 
number of terms which grows with L. Much of our effort has been di- 
rected toward finding good, easily computed bounds and asymptotic 
formulas. 

The probabilities of having exactly A^ coincidences are also obtainable 
but they have more complicated formulas. A detailed derivation is given 
only in Part II. 

In Part IV, we consider the probability of no coincidence in higher 
dimensional problems. The methods of Parts I-III fail in higher di- 
mensions, but we are still able to derive some bounds. An exact formula 
is derived for the probability of no coincidences within a single two- 
dimensional Poisson pattern in a rectangle with sides ^ 25. We also give 
particular attention to coincidences in a three-dimensional cylinder. 

Part V contains numerical results. 

Heduction of the Exam-pies to the Theory 

We now msh to see how answers bearing on the practical problems 
previously listed may be found from this study. 

The literature on Geiger counters (see bibhography in Feller^) is con- 
cerned with statisticH of the number of counts registered In a given long 
time, t. The basic problem is to test the hypothesis that the particles 
arrive in a Poisson sequence. To this problem, then, are relevant the 
formulas for the probability of A'^ coincidences in one pattern given in 
Part II, and the bounds and asymptotic results there derived. 

The problem of coincident flaws in an electric cable is three-dimen- 
sional, and we have various approaches leading to the probability of no 
coincidences which are vahd under different circumstances. If the cable 
contains only two uires (with possibly different flaw densities), then the 
problem reduces to the one-dimensional (^asc of coincidences between 
two Poisson patterns treated in Part I. If the diameter of the cable is 
small with respect to 5, and if the density of flaws is the same on each of 
the n wires in the cable, we have the situation of n identical patterns 
treated in Part III. If, in addition, n is very large, we may ignore the 
fact that coincident flaws on a single wire do not cause short circuits, 
and think of coincidences \\ithin a .single pattern (Part II). Without the 
assumption that the diameter of the cable is small with respect to 6, the 
problem is no longer reducible to a one-dimensional foim. Section 4.4 
Is especially devoted to tliick cable, and to ])roducing a lower bound for 
the probability of no coincidences in this three-dimensional situation. 

The literature on Poisson patterns in a line segment contains the fol- 



1008 THE BELL SYSTEM TECHNICAL JOURNAL, JULY 1957 

lowing related papers. C. Domb' finds the distribution function for the 
total length of the set of points lying within distance 5 of a pattern point. 
P. Eggleton and W. O. Kermack'' and also L. Silberstem^ consider ag- 
gregates, which are sets of k pattern points all contained in an interval 
of length 5. In the special case k = 2, aggregates are our coincidences. 
These authors find the expected number of aggregates but not the prob- 
ability of A'" aggregates. 

1 CONCIDENCES BETWEEN TWO PATTERNS 

1.1 Integral Equation 

Consider two Poisson patterns of points on the real line, the first with 
density X (points per unit length) and the second with density /i- We want 
the probability F{L) that in the segment from to L there is no coinci- 
dence between a point of pattern No. 1 and a point of Pattern No. 2. 
F{L) will be formulated in terms of the conditional probabilities 

Pi(L) = Prob (no coincidence, given Pattern No. 1 has point at L), 

PiiL) = Prob (no coincidence, given Pattern No. 2 has point at L). 

If L ^ S, Pi{L) and PiiL) are the probabifities that patterns No. 2 
and No. 1 are empty: 

Pi(L) = e-''^ P-tiL) = e^'\ if L ^ 5. (1-1) 

If L > 5 and Pattern No. 1 contahis a point at L, there are two ways 
that no coincidences can occur. First, Pattern No. 2 may fail to have any 
points anywhere in the interval [0, L\. The probability of this event is 
exp - fxL. The second possibility is illustrated in Fig. 1 (using circles for 
points of Pattern No. 1 and crosses for points in Pattern No. 2). Pattern 
No. 2 has points in (0, L)\ the one closest to L is at y < L - 5. Since 
the interval {y, L) contain.s no points of Pattern No. 2, the probability 
of finding this closest point, //, in an interval, dy, is 

exp l-fiiL - y)h dy. 

The interval {y, y + S) must be free from points of Pattern No. 1 (prob- 

NO COINCIDENCES EMPTY NO CROSSES 
^ . 1 .^ -J ^^— ' — ^ 



■^(r^ @ ^« 1 ^t^^ ® 



Fig. 1 — Patterns without coincidence. 
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ability exp (~U)) and the interval (0, y) must contain no coincidences 
{probability P^iy))- One obtains finally 



PM = e 
and similarly 

P-2{L) = e 



.-Ml- 



l+/.e-^^jf' \'"'PMdyj, (1-2) 



..L-S 



l+Xe-"^!' ' e^^P.iy) dy'j. (l-'S) 



The solutions Pi{L) and P.iL) are determined uniquely by (1-1), 
(1-2) and (1-3). For (1-1) determines them for g A ^ 5 and the in- 
tegrations indicated in (1-2) and (1-3) will pro\-ide the solutions in g 
L ^ (n -\- 1)6 when they are known in ^ Z. ^ t^S. Pi(L) and PziL) 
are piecemse analytic; the analytic form of the solution changes 
each time L passes an integer multiple of 5. These analytic expressions 
soon become complicated and are less useful than the bounds and ap- 
proximations given later on. 

To compute F(L), consider the last place l^efore L at which either Pat- 
tern No. 1 or No. 2 has a point. The probability that this last pomt lies 
between x and .i; + dx and belongs to Pattern No. 1 is exp [ — (X + fi) 
(L - x)]\ dx (Fig. 2). This term multiphed by Pi{x) and integrated from 
to L gives the probability of no coincidences if the last point is a circle. 
A smiilar integral gives the probability if the last point is a cross. Finally 
there is probability exp [-(X -\- y.)L] that neither pattern has a last 
point [i.e., (0, L) empty]. Then 

F{L) = e-'^+^'^ \j + I' e*^+'"^(XA(.r) + t.P,{x)) d^. (1-4) 

1.2 Solution by Laplace Transforms 
For t" = 1 or 2, let 

Piis) = [ PiiDe-"- dL. (1-5) 

Jo 

Replacing PiiL) in (1-5) by (1-1) for ^ L ^ 5, by (1-2) for 5 g L, 

NO COINCIDENCES EMPTY 



-^^ ^ '•(r^ ro i O i 

X x + dx 
Fig. 2 — Patterns without coincidence. 
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and interchanging the order of integration of a double integral, 



Similarly, 



so that 



and 



(s + \)p2is) = 1 + Xe^"+''"'"'Vi(s), 



, , s + A + tie ^^ /, px 

^'^'^ = (. + X)(s + m) - Xm^-^'^^"^"" 



/ -, _ s -\- fi -\- \e Q_y-, 

P'-^^> - (s + x)(s + m) - XMe-^'^+'^'»* 

Likewise, using (1-4), the Laplace transform /(s) of F{L) is 

.. -, _ 1 + Xpi(s) + Ilp2{s) 
^^'^ ' X + M + S • 

As one might expect from the piecewise analytic character of Pi(L) 
and P2{L) there is no convenient way of transforming /(s) back to F(L). 
By evaluating residues of /(s) exp (sL) at the poles of /(s) one might ex- 
press F{L) as an infinite series of exponential terms. The most slowly 
damped term in this series can be expected to approximate F(L) when L 
is large. The poles of f{s) are at the zeros of the denominator D{s) of 
pi(s) and p2is): 

D{s) ^ (s + \){s + m) - \ne-''''^''^"\ {1-9) 

Since D{x) > for a: ^ and both D(-\) and i)(-M) are negative, 
it follows that /)(s) has a real zero s = —a with a < Min (X, m)- 

The zero s - -a of D{s) is the one with the largest real part. For, 
letting s = a; + iy, we have in the half plane x '^ —a 

I (5 + X)(6- + m) I - AM I e I 

= I s + X I ■ I s + M I - \p.e-"-'^^'^^'' 

^ (.^ + x)-(x + m) - XMeT'*'^"-^''* ^ 0. 

Also, if // 7^ the ^ sign in the above proof cun be replaced by > and 
one concludes that all other zeros of D{s) = satisfy 

Res < -h 
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for some b > a (note that the left hand side of the preceding inequaUty 
does not approach asy approaches ± =o ) . 

The pole of /(s) at s = —a contributes to F(L) a dominant term 

F{L) ^ X^ + M^ - (X + /^)a + 2Xm^"^^-^--°" .^ jq. 

^ ^ (X -I- M - a)[X + M - 2o + 25(X - a)(M - a)] ' ^ ^"^ 

In (1-10) the error is 0(exp — bL) for large L. 

When 5 is small, we find a = 2Xm6 + 0(5^) and (1-10) becomes 

F{L) ;^ [1 + 0(5')] exp - [2X^5 + 0{5')]L. (1-11) 

It is interesting to note that a simple heuristic argument also leads to a 
formula like (l-U). When 5 is small and L is large, one expects that the 
intervals of length 28 which contain points of Pattern No. 1 at their cen- 
ters will comprise a total length near (XL) (26) of the line segment (0, L) . 
The probability that a set of length 2XL5 shall be free of points of Pat- 
tern No. 2 is exp — 2X/i5L. 

1.3 Bounds 

In this section we derive some relatively simple expressions which are 
good upper and lower bounds on F{L). Both bounds have the same func- 
tional form: 

K{A, B- L) = 24±i^ ,-- + (i _ X^_±i^) ,-(^+M)^ (1,12) 
X + M — « \ X+ju— a/ 

In (1-12), a is again the smallest real solution of D( — a) = 0. ^ and B 
are positive constants which are related by 

'^ ^ M ^-(\+M-o)i ^ X — g ^(x+p-a)s d-lSl 

B fi — a X 

K{A, B; L) becomes an upper bound or a lower bound depending on ad- 
ditional restrictions which will i>e placed on .4 and B. 

To get the lower bound, we restrict .4 and B by the inequalities 

A < e'"-'^', B < f'"-''^ (1-14) 

and 

A <(l-^c^\ B<(l-^^c'r (1-15) 

We first prove that (1-13), (1-14), and (1-15) imply 

■ Pi(L) > Ae-"'', P,{L) > Be^"''. (1-lG) 
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When ^ L ^ 5, (1-16) holds because of (1-1), (1-14), and the in- 
equaKties a < \, a < fi. If (1-16) were not true for all L there would be 
a smallest value, say L = X > 8, at which at least one of the inequali- 
ties (1-16) would become an equality. Suppose the inequality (1-16) on 
Pi(Z) fails. Using (1-16) for L < X, and (1-2), 

/ Jp-aHX-S) ■, 



Pi(X) > e-"^ 1 + S^e 

\ n - a 

1 _B^^^)e-- by (1-13), 

> Ae'""" by (1-15). 

This contradicts our assumption that (1-16) fails for Pi{X). A similar 
proof shows (1-16) cannot fail for P2iX). 

Having proved (1-16) we now substitute these bounds into (1-4) and 
integrate to get F{L) > K{A, B; L). 

To make (1-12) into an upper bound it is only necessary to replace 
(1-14) and (1-15) by 

^ > 1, B> \, (1-17) 

and 

A>(l-^ e"', B > (l - ^) ^\ (1-18) 

The proof that now F(L) < K{A, B; L) proceeds exactly as before 
but with all the inequality signs reversed. 

Both bounds are dominated by an exponential term exp — aL, as is 
the asymptotically correct formula (1-10). In typical numerical cases the 
coefficients multiplying this term in the three formulas agree closely. A 
numerical case is given in Part V. 

1.4 Probahilitij of N Coincidences 

The methods of Sections 1.1 and 1.2 can also be used to find the prob- 
ability Fn{L) that there be exactly A'' coincidences in the interval (0, L). 
It might appear most natural to define A'^ to be the number of pairs of 
points (.r, z), x from Pattern No. 1, z from Pattern No. 2, such that 

\x - z\ < 8. (i) 

However, we add the additional requirement that x and z be "adjacent" 
points; i.e. 

the interval (x, z) is empty. (ii) 
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For example, in Fig. 8, we would count A'^ = G coincidences even 
though there are 18 pairs which satisfy (i). In cable problems it appears 
reasonable to count coincidences as above. If we assume that all flaws 
are equally bad, then a short circuit is likely to develop only across an 
adjacent coincidence; our A^ is the number of place.s on the cable at 
which a short circuit can form. Another interpretation is that the cable 
can be cut into exactly iV + 1 pieces each of which contain no coinci- 
dences. 

Let Pi,n{L) be the conditional probability of having A'' coincidences 
in (0, L) knowing that there is a point of Pattern No. 1 at L. The Lap- 
lace transform of Pi_n(L) turns out to be the coefficient of t'^ in a generat- 
ing function of the form 



(X + s)(m + s) -XmS^-' 
where Q = e"*^'''''+''' (i _ /j _|_ /^ Interchanging X and m one gets the gen- 
erating function psit, s) for the Laplace transform of the probability 
P2,n(L) of N coincidences, given a point of Pattern No. 2 at L. Finally 
the Laplace transform of Fn(L) is the coefficient of t' in the generating 
function 



fit, s) = 



X + A» + S 



Since /(^, s) is a rational function of (, it is easy to find the coefficient of 
l^. The poles of this fmuition are again just zeros of D{s). Now, however, 
the poles are higher order poles. For large L an a.symptotic formula for 
Fn{L) has the form exp — aL times a polynomial in L with degree de- 
pending ou N. 

For more details about this method we refer the reader to Part II 
where a similar, but less involved, calculation is carefully done. 

II SELF-COINCIDENCEg IN ONE POISSON PATTEHX 

^.1 Integral Equation 

In this part we shall consider a .single one-dimensional Poisson pattern 
with density X and ask for the probability i^,v(/.) that in the interval 
(0, L) the pattern have exactly N coincidences. We count coincidences 

I ^Ttfc ^ 1 @^^ K X • > *- ' 3 «-@ 1 

1^ L 

Fig. 3 — Patterns with six comcidenees. 
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as in Section 1.4; a pair (x, z) of pattern points contributes one coinci- 
dence to the total number A^ only if both \x — z\ < 5 and the interval 
between x and z is empty. 

Note that Fq{L) is related to the distribution function for the mini- 
mum distance between the points of the pattern in (0, L) : 

Prob (min. dist. ^ a) = 1 - Fo(-L), 

where it must be remembered that /''o(i) is a function of 5. 

As in Part I, we first define the conditional probabilities Pn{L) = 
Prob (exactly A^ coincidences in (0, L), given a point at L). We then 
have the follo^-ing equations: 

liL^b, PsiL) = ^ e-'\ all N. (2-1) 

liL^6, FAL)=^^^_e-'\ N^X. (2-2) 

If L > 5, and N ^ 1, the probabihty of exactly N coincidences in (0, L) 
equals the probability of N coincidences up to the last point of the pat- 
tern in the interval (0, L) - and if there are to be any coincidences, there 
must be points of the pattern in (0, L). Hence, if L > 5, iV ^ 1, 

FAD = C PAL - i/)e-'''X dij. (2-3) 

Jo 

li N = 0, the same argument applies, but there is also the possibiUty 
that there are no points at all of the pattern in (0, L). Hence, if L > 5, 

F,{L) = e-"' + f ' PoiL - y)e-'"\ dy. (2-4) 

Now let us consider the case where there is a point of the pattern at L. 
Then if the last pomt preceding L is between L - 5 and L, this point 
and the point at L will create a coincidence; if there is no point within 
(L - 6, L), then all coincidences are within (0, L - 5). Hence, if L > 5, 
andiV ^ 1, 



PAh) = f P^^iiL - y)}.e'"' dij + e-''FAL - 8). 

Jo 



(2-5) 



For the case A'' = 0, we cannot allow a point in the interval {L - 8,L), 
and hence, if L > 5, 

P,(L) = e-'^'F.iL - 5). (2-6) 
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^ J Laplace Transform of Fn{L) 

To analyze the system of equations which is given by relations (2-1) 
through (2-6), we introduce the generating functions 

f{L,i) = ZFAQt'', 

■v-o 

and 

p{L,i) = J:PAL)tr 

If L > 5, we obtain from (2-3) and (2-4) the relation 

e'V(L, - 1 + / p{w, rte'^X dw, (2-7) 

Jo 

and from (2-5) and (2-6) the relation (again if L > 5) 

e^'-piL, I) = U [ p(w, Oe^'" dw + e'"'"*'/(L - 5, 0- (2-8) 
Jt-i 



If we differentiate (2-7) and (2-8) with respect to L, and then apply 
(2-7) differentiated to simplify the last terms of (2-8) differentiated, we 
obtain, still only for L > S, 

f'{L, t) + X/(L, ^ \p{L, I), (2-9) 

p'(L, t) + \(1 - i)p{L, t) - \e~^\\ - ()piL - B, t). (2-10) 

It is easy to check from (2-1) and (2-2) that if Z/ ^ 5, then 

/T i\ — xm— 1) 
p{L, {) = e 

and 



;(L,0=e-"((^^ + l), 



and hence (2-9) is vahd for all L, but the left side of (2-10) vanishes if 
L ^ 6, Hence we may take Laplace transforms of (2-9) and (2-10). If 
we define 



and 



A(s,t) = ( f{L,t)e-''dL, 
Jo 

Ii{s, t) = [ p{L, l)e-'" ilL, 
Jo 
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we obtain from (2-9) , which we now know to he valid for all L, 

(X + s)A{s,t) - 1 = \BisJ), (2-11) 

and from (2-10), by recalling that the left side vanishes for L ^ S, 

sB{s, - 1 + X(l - l)B(s, t) = X(l - Oe"^''"'"*B(s, 0- (2-12) 
Hence 

^fa '> - . + x(i - o[i - .^'-^"r ^'-''^ 



and 



A(s, t) = -J— (1 + XB(s, 0). 

A + S 



If we denote the Laplace transforms of Pn(L) and Fn{L') by pf,(s) and 
/«(s) respectively, then 



and 



Us) = -^ (Xpo(s) + 1), 
X -p s 

Uis) ^^^pAs) for iV = 1,2, 

A -p S 



(2-15) 



2.3 Exact Formula for Fd{L) 

It is possible to solve (2-1) through (2-6) in piece\vise analytic form by 
computing recursively from each interval of length 6 to the next one. We 
shall obtain the piecewise analytic form for Fa{L) by a direct derivation 
essentially due to E. C. Molina.* 

Suppose k is the number of pattern points which fall into (0, L). Let 
Xi denote the distance between the i - 1^* point and the i* point (xi is 
the distance from to the first point) as shown in Fig. 4. The configm-a- 



Xi. Xk 

L . - L. 



-@-@ ■^±rrti ^^ 't^ @ @- 

12 1-1 I ^\ 

Fig. 4 ^ Definition of Xi . 
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tion of points 1, • • ■ , k on the line is represented by a single point 
{xi, ' • ■ yXk)\n the polyhedron T in fc-dimensional space defined by the 
inequalities 

T: g a-i , • ■ • , ^ atfc , xi + x^ + ■ ■ ■ ^ Xk ^ L, 

and the probability distribution of the point (a-j , ■ • ■ , xC) in T is uni- 
form. The configurations with no coincidences lie in a smaller polyhedron 
T' consisting of all points of T for which 5 ^ .T2 , ■ ■ ■ , 5 ^ .t*, . Given k, 
the conditional probability that there be no coincidences is the ratio of 
two /i>dimensional volumes Vol (7")/VoI {T). 

Vol (TO = if L g (A: - 1) 5. 

For larger values of L let //i ^ Xi , 1J2 ^ X2 — S, y^ = xx —8, • • • , Uk = 
Xk — 6. Then 7" becomes a polyhedron of the form 

T": g 7/1 ,0 ^ //, , ■ ■ -,0 ^ ?/, , l/i-\-y2--- +yi-^L~{k- 1)5. 

Since the transformation from .t's to y's has determinant equal to one, 
T" has the same volume as T'. However, T" is now seen to be similar 
to T but mth sides of length L — {k — 1)5 instead of L. The volume 
ratio sought must be 



/ L-jk- Ds V 



Since k has the Poisson distribution Tvith mean \L we obtain finally 

The piecewise-analytic character of Fo(L) is eWdent; increasing L by 
an amount 6 increases the upper limit on the sum by one and thereby 
adds a new term to the analytic expression for F{L). 

2.4 Asymptotic Formula for Fx(L) 

Similar exact formulas could be found for all the Fs{L), but they are 
both comphcated and inconvenient for computing if L/5 becomes large. 
It is thiis natural to aim for asymptotic results and for bounds connected 
with them. 

The Laplace transform of Fy{L) is given through (2-14) and (2-15) 
above. The pole of fy{s) with largest real part is a pole of order A'^ + 1 
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at a real negative point 

s = —a > —\. 
For large L, the asymptotic behavior is given by 

tsKi^) - ^^ _ ^^^^ ^ g^^ _ ^^]^j |_j_ ^ 5^^ _ „jj . 

where the error term is 0(L""^ e"°^) if A?" ^ 1. Such a formula, then, is 
a good approximation for fked N as L increases; for fixed L, however, 
it will fail to be good for sufficiently large A''. 

Ji N = 0, the asymptotic form is 

^'^^^ ^ (X-a)[l + 5(X-a)l'~" 

but the error term now decreases at a more rapid rate, as may be seen 
by including the contributions of some of the complex poles of /o(s). To 
find these poles, set 

s + X = Xe-''+'**. 

If 

s = — X + r exp (id), 

one obtains the simultaneous real system 

2wm. - e = 6r sin 6 (m integer), 

log (r/X) - -5r cos 6. 

The first equation defines an infinite family of curves in the s-plane (see 
Fig. 5). The second equation defines a single curve which intersects the 
family at poles of p{s). 

2.6 Bounds on Fa(L) 

As in Part I, we may derive bounds on Fo(L) from the integral equa- 
tion, and obtain 

Since a = X"5 + 0(5^) for small 5, the bounds are very close if \8 is not 
too large. 
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Fig. 5 — Solution of s + X = \c 
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III COINCIDENCES BETWEEN 71 POISSON PATTERNS 

3.1 Integral Equation 

In this part we consider n one-dimensional Poisson patterns and ask 
for the probability, F{L), that in the inter\'al (0, L) no pair of points 
from different patterns are coincident. Unlike Part I, we now consider 
only the case in which aU n patterns have the same density X. Let P{L) 
be the conditional probability, given that Pattern No. 1 has a point at 
L, that there are no coincidences in (0, L). 

If g L ^ 5, P{L) = exp - (n - \)\L, 

If 6 < L, 

P{L) - ,-<"-^'^^ U^ {n- l)\e-^' f e^"'''^"P{y) dij\ 

by the same sort of argument used in Part I. Then F(L) will be given 

by 
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S.2 Bounds and Asymptotic Formula 

The Laplace transform of P{L) is 

p(s) = {s+ (n - 1)\(1 - c^'"^-^^^^)}-' (3-1) 

which has one real pole at a negative point s = —a, a < {n — 1)X. 
Again it is this pole which contributes the dominant term to both P{L) 
and F{L) for large L. We find 



F(L) 



n\e 



(1 + [{n - 1)X - a]5)(nX - a)' 
To bound P(L) by expressions of the form A exp( — aL) one finds 
that A > 1 wiW give an upper bound and 

will give a lower bound. The corresponding bounds on FiL) are of the 
form 

n\ — a) n\ ~ a 

3.3 Exact Solution 

As in Part II an exact formula for F{L) may be given as a finite sum. 
We now derive it from the Laplace transform, 

/(s) = {s + nKY' (1 + nkvis)), 

of FiL). We may use (3-1) to expand /(s) into the series 

The identity 

(s + n\)-\s + {n - DX)-'-' 

1 -^ , . , (3-3) 

provides a partial fraction expansion for the k''^ term of the series (3-2). 
Transforming (3-2) term by term with the help of (3-3) we find 

F{L) = e-"'^[-(n - 1)]'^'*'+^ 

This is the desired formula for F{L). 
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IV MULTIDIMENSIONAL PROBLEMS 

4.1 Two-Pattern Lower Bound 

We now derive some results on the probabilities of no coincidences in 
some multi-dimensional situations. The simplest one is a lower bound 
for the case of two Poisson patterns. 

Theorem: Consider a d -dimensional region of volume V containing two 
Poisson patterns ivilh densities \ and fi. Let S{8) be the volume of the d-di- 
mensional sphere of radius 5. The probability of no coincidences between 
the two patterns has the lower hound 

Proof 

Let the pattern vnih density X be called the X-pattern and the other 
the /i-pattern. Given any X-pattern of /,- points there will be no coinci- 
dences provided only that a certain region 2' contains no points of the 
jLt-pattern. T consists of all points of the ^-olume V which he in any of 
the spheres of radius 5 ceTitered on the A- points of the X-pattern. Since 
these spheres may overlap and may extend partly outside the volume 
V, we have 

volume oiT S k 8(5), 

and 

Prob (no coinc, given /.■ points) = exp ( — ^ volume of T) 

^ exp(- l-^S{&)). 

Since the number, k, of points of the X-pattern has the Poisson distribu- 
tion with mean XT'" the (unconditional) probability of no coincidences 
has the lower bound 

E(xy) -w -kiis(.s) 

Summing the series one proves the theorem. Interchanging X and ^ in 
the theorem gives another lower bound. The one stated above is the 
better of the two if X < n. 

The difference between the lower bomid and the true probability 
comes from two sources: (a) The overlap between the k .spheres; this 
will be a small effect if \'S{2S)V is small, and (b) the spheres which 
extend partly outside the volume V; there will be relatively few such 
spheres if only a small fraction of the volume V lies within distance 5 
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of its boundary. Hence in some cases the lower bound will be a good 
approximation to the correct value. 

It may also be noted that no real use was made of the spherical shape 
of the volumes S{S). If one wants to consider a point of the M-pattern 
to be coincident with a point of the X-pattern if it lies in some other 
neighborhood, not of spherical shape, the same lower bound applies 
but with S{5) replaced by the volume of the neighborhood. 

4.S Single-Pattern Lower Bound 

A similar derivation in the case of a single Poisson pattern leads to: 

Theorem: Ld a Poisson pattern of density X be distributed over a d-di- 

mensional region of volume V. Let 8(6) be the volume of the d-dimensional 

sphere of radius 5. Then the probability of no coincidences is at least as 

large as 

- e-'^'U + X.S(5) }""'". 

The theorem will follow from another bound which is slightly more 
accurate but much more cumbersome. 

Lemma 

In the above theorem a lower bound is 



~-\v 
e 



( 1 + XV + E ^ n [1 - js{B)/y] ■ (4-1) 

\ t=2 A'! 3 = 1 / 



Proof of Lemma 

The probability sought is of the form 

E "-"' ^ P^ , (4-2) 

where pk is the probability that, when exactly A- points are distributed 
at random over V, there are no coincidences. To estimate pi- , imagine 
the k points to he numbered 1, 2, ■■■/(■ and placed in the region one at 
a time. If no coincidences have been created among points 1, ■■■, j 
(which is an event of probability pj) the probability that the addition 
of point j + 1 creates no coincidence is just the probability that this 
new point lies in none of the j spheres of radius 5 centered on points 
!, ■■-,/. The union of these j spheres intersected with the volume T' 
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is always of volume ^ jS{6). Hence 

P.>. S pAI - JS{8)/V], 
or 

p. ^ n [1 - jSi5)/V]. (4-3) 

When (k — l)S{d) > V the above argument fails because the later 
terms of the product are negative ; in this case we use the trivial bound 
Pk ^ 0. Combining (4-2) with (4-3) the lemma follows. 

Once more the bound may be expected to be almost correct if }^VS{2S) 
is small and if most of the region V lies farther than S away from its 
boundary. The bound is also correct for non-spherical neighborhoods 
(see discussion of previous theorem). 

When V/S{5) is large, the sum (4-1) is unwieldy. If we let H equal 
V/Si8), we maj' rewrite the typical term in the sum as 

^n (1 - J/H) - ^^^^^H(H -1)...{H- k + 1). 

If H happens to be an integer, this equals 



^. , i\V/H)\ 
so that the complete sum (4-1) equals 

We will no^^' prove that if H is not an integer, the sum always exceeds 
(4-4), so that (4-4) is a lower bound in all cases. We wish to prove that 

IH]+1 k 

1 + E TT.HiH - 1) ...(//- A: + 1) ^ (1 + aO" (4-5) 

t=l K\ 

for any positive H, in which event the theorem follows mth 

X = ^ and H = V/S{5). 

The inequality (4-5) will be proved by induction on \H]. If \H\ = 0, 
then we are required to show that 



1 + Hx ^ (1 + x)" 
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for ^ i/ < 1. This follows immediately from the concavity of 

Suppose now that (4-5) holds for a value H. If we mtegrate both 
sides of (4-5) from to x, we obtain 

which may be rewiitten as 

1 + E n (^ + l)t^) ■ ■ ■ (H - /c + 2) ^ (1 + .r)"+\ 

(.-=1 'i^! 

This completes the induction, and the proof of the theorem. 

4.5 Another Lower Boimd {Any Number of Patients) 

Another kind of lower bound can be derived which sometimes will 
be better than the above bounds when the region V has a large fraction 
of its volume mthin 6 of the boundaiy. For example, V might be a 
three-dimensional circular cylinder (a cable) with a I'adius which is com- 
parable to 5. 

To derive this bound one first finds the expected number, E, of co- 
incidences in V. An upper bound on E will also suffice. Then it is noted 
that 1 — £^ is a lower bound on the probability of no coincidences. For 
if Qif is the probabiUty of finding A'' coincidences, 



AT-l 



44 Thick Cable 

For example, we now give a lower bound which is of interest in con- 
nection with the problem of a cable with many wires. 

Theorem: Let a Poisson pattern of points with density X be plmed in a 
cylinder of length L and radius R > 5. The probability of finding no co- 
incidences in the cylinder is at least as great as 

Proof 

Introduce cylindrical coordinates r, <p, Z so that the cylinder is de- 
scribed by 

r ^ R, ^ Z ^ L. 
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Consider first any pattern point (r, p, Z) with Z-coordinate satisfying 
S ^ Z ^ L — 5. Let arrows be drawn from this point to all other pat- 
tern points (if any) within distance 5. The expected number of arrows 
draAMi from this point will be \G{r) where G(r) is the volume of the 
intersection of the cylinder mth a sphere of radius 6 centered at the 
point. For points near the ends of the cylinder (Z ^ 5 or L — S ^ Z), 
the expet^ted number of arrows will be less than \G{7-). Since the proba- 
bility of finding a pattern point in a httle volume element dV is XrfV, 
wo conclude that the expected number of arrows drawn in the entire 
cylinder will be less than 



/// ^"«fr) "''■ 



cylinder 



If the cylinder has A'' coincidences, there will be 2iV arrows (each point 
of a coincident pair appears once at the head of an arrow and once at 
the tail). Hence the expected number of coincidences is 



E 



^ xVL [ G{r) r dr. (4-7) 

Ja 



Since an exact formula for G{r) is rather cumbersome, we are content 
with a simple but close upper bound. Jl r ^ R — S then clearly G{r) = 
47r5'/3. If r > ff — 5 we get an upper bound on G(r) by computing the 
shaded volume in Fig. 6; the intersection of the sphere with a half -space. 

Gir) S [25' + 3(K - r)5' - {R - rf] vr/S. 

Substituting these expressions for G{r) in (4-7), integrating, and using 
(4-6) the theorem follows. 

The approximation to G(r) whidi was made above is bad when R is 
much less than 5, but in this case good estimates may be obtained from 
the one-dimensional results of Part II. Note also that if X is large enough, 
the bound becomes negative and is therefore useless. 



-SPHERE 




Fig. 6 — A region for estiniiitiiig G(r). 
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4.6 Upper Bounds 

Good upper bounds appear even harder to get than lower bounds. 
One procedure is to divide the region V into a number of smaller cells. 
If each cell has probability, p, of no coincidences and if there are K 
cells, then p^ is the probability of no coincidence in any cell. If there is 
no coincidence in V there will be none in any cell; hence p^ is an upper 
bound on the probability of no coincidence in V. 

Of course, p^ is too large because of the possibility of a coincidence 
between two points in different cells. It follows that p^ will be a close 
bound only if the cell size is made large; but then p becomes hard to 
compute. 

For example, consider self-coincidences in a single Poisson pattern in 
a large region of area V in the plane. Cover this area with an array of 
hexagonal cells of side 5/2 as shown in Fig. 7. The area of each hexagon 
is 3 V3 ^VS so the number of cells used mil be about K = 87/3\/3 S^. 
A cell has no coincidence if it contains at most one pattern point, hence 

p = (1 + X8\/3 SVS) exp - 3-s/3 XS'/S- 



The upper bound is 



K -w 

p = e 



(i + ^ ^^,y.™v3« 



which has an interesting resemblance to the lower bound 

e (1 + TTAd ) 




Fig. 7 — Pattern for studying coincidences in a plane region. 
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4-6 An Exact Calculaiimi 

The upper and lower bounds in Section 4.5 are not very close, largely 
because of the small size of the hexagonal cells. An improved upper 
bound may be obtained using square cells of side 25. We can calculate 
■p for small rectangular cells but only if we redefine our notion of coin- 
cidence in terms of square neighborhoods instead of circular neighbor- 
hoods. That is, points {xi , iji) and (.ra , ^2) are now considered coincident 
if simultaneously 

I -ri - .1-2 1 g 5, and \yi — yi\ ^ S. 

The result we get is the only exact calculation of a non-trivial multi- 
dimensional coincidence probability knouai to us. 

Consider the rectangle ^ x ^ L, ^ y ^ M with L and M both ^ 
25. If L is less than B, two points are coincident if and only if their y-co- 
ordinates dilTer by less than 5. The problem then reduces to a one-di- 
mensional coincidence computation such as we ga\'e in Part II. There- 
fore, suppose both L and M are greater than 5. 

There is probabiUty 

(XLM)' ^..LM 

that the rectangle contams k points. We therefore subdivide the problem 
into cases of the form "given /,-, find the probability that the k points 
have no coincidences". Only five of these cases have a non-zero answer. 
To show this, divide the rectangle into four rectangles of sides L/2, M/2; 
if A; ^ 5 one of the.se rectangles must contain more than one point, and 
so a coincidence. The remaining cases k = 0, 1, 2, 3, 4 may be further 
subdivided according to which pairs of .r-coordinates are less than 5 
apart. Let us number the k points (.I'l , yO, ■ ■ ■, (.rt , yk) in such a way 
that the ;r-coordinates are in order Xi ^ .r2 ^ ■ ■ ■ ^ Xk . If, for some i, 
Xi+2 ^ a:, -\- 5, then the subcase in question contributes zero to the 
probability of no coincidences because all of | x; — .t,-+i |, | a.',-+i — .r,>2 |, 
I .r,+2 - Xi I are ^ b and at least one of | //, - ^,+1 [, [ ?/,+! - //,+2 |, 
I yi+2 — Vi I is ^5. The only subcases which remain to give a non-zero 
contribution are the nine Hsted in Table I. The number in the "subcase" 
column is k. The uext column contains the a:-inequahties which define 
the subcase. The probability that the k ordered .^-coordinates satisfy the 
stated inequalities is listed as probj . If the .r-inequalities are satisfied 
there will be no coincidences if and only if | //& — ^a | > 5 for every m- 
equality \xb — Xa\ ^5 given in the .r-inequahty column. These y-m- 
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Table I 



Subcase 


X inequ. 


probi 


y inequ. 


prob]/ 





— 


1 




— 


1 


1 


— 


1 




— 


1 


2(a) 
2(b) 


X2 — Xi > S 

xi - Xi ^ S 


(1 - 5/Ly 

25L - 5^ 




1 ya - yi 1 > 5 


1 

(1 - s/My 


3(a) 


X2 — Xl ^ S 

X3 - Xi > S 


1(1 - 6/Ly 




1 2/2 - yi 1 > 5 


(1 - 5/ My 


3(b) 


X2 — X| > 5 


j(i - 5/Ly 




1 ya - ys 1 > 3 


(1 - s/My 


3(c) 


X2 — ii S 5 
X3 - X- ^ d 
Xj - Xi > 5 


j(l - &/Ly (^^ - 





1 ya - ya 1 > 5 
1 yi - yo I > 5 


Ui - s/My 


4(a) 


l3 - 12 ^ 5 
X] — Xi > 5 
Xi — Xi > S 


Ml - s/Ly 




1 ya - yi 1 > 5 
1 ys — y2 1 > 5 
1 !/* - y^ 1 > s 


■{•,a - 5/MY 


4(b) 


X3 - X2 > S 


Ki - 5/Ly 




1 y2 - yi 1 > 5 
I y4 -- y3 1 > 5 


(1 - 5/''^)' 



equalities are listed in the third column and the prol^abilitics that they 
are satisfied are listed as prob„ . The probability of no coincidences is 

X^ Qk probi proby 
where the sum is over all nine subcases. The sum is 

exp {-\LM) jl + \LM + |- {LV- - 5'(2L - 5)(2ilf - 5)] 

+ ?^ (L - 5)-{M - b)\2LM + L5 + MS - 45') 

at 



+ ga-a)<(M-«j. 



If L = M = 25, this reduces to 
exp (--IS'X) 



l + 45^+^5V+i^5V + 5^ 



^ + ^4^^^^]- 
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A sample of one of the above computations may be instructive. Con- 
sider, for example, Case 4(a). We have ^ .ri ^ x-i ^ .7:3 ^ .u ^ L, 
and require: 

;l-3 - .1:2 ^ 5, 
-Ta - Xi> 5, 

.n - a-2 > 5. 
The probability of this is 

/L [-l.— i nTJ— i ~L 

/ / / dxi dxi dx2 dxz 

_3=J •'j-j^ijra— J •'j;i=0 ''^4=12+! 

- Yll \ {L - 1-2 - b){xz - 5) dx^ dx2 



8 (L-sr^l/ _^Y 



L* 24 3\ L 

In the ^/-direction we requii'e \ 1/2 — yi\ > 5, \ y^ — y^l > 5, \yi — 
//a I > 5, and there are no order restrictions. Assume first that y-i < y-.. . 
Then the probability that yi and y^ satisfy their restrictions is 

/ j/3 - h \ ( M - yn- 8 
\ M )\ M 

Hence, the probability for satisfying all the conditions is 

/•■" r'" ( y^ - A Av -y.-b \ dihdm^_^(. _±S 

J, Jo \ M J\ M /MM 24 \ Mj ' 

Interchanging yz with ys and yi ^'ith */4 shows that the assumption 7/2 > 
//3 yields the same answer, so that the required probability is 

12 V M, 



V NUMEIilCAL WORK 

5.1 Coincidences between Two Patterns 

5.1.1 Machine Com-piitation of F{L) 

To compute the probability of no coincidences in a line of length L 
directly, it is convenient to transform equations (1-2) through (1-4) into 
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the following differential difference equations: 

fO if .1: g S 
P2(.r - 5)fxe 



if .T > 8, 



Pa'Cr) + XP2W = 



[0 if :r ^ 5 



[P,{x - 5)X.-"'^'"* if x > 8, 
F'ix) + (X + n)F(x) - XPiCr) + fiPiix), 
p,(0) = p,(0) = P(0) = 1. 

These have been solved on a general purpose analog computer mth 
the aid of a lumped-element approximate delay line for a number of 
cases. We have chosen for illustrative purposes the parameters X = 5, 
fi = IQ, 8 = 0.02, and L g 1. The exact solution, together with various 
approximations to be described in the sequel, is plotted in Fig. 8, where 
the exact solution is labelled yi . 
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Frc. 8 — Probability of no coincidences between two one-dimensional Poisson 
patterns with \ = 5, /i = 10, if 5 = 0.02. 
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5.1.2 The Asijmptotic Foi-rmila 

All approximation to the probability F{x) of no coincidences is given 
by the asymptotic formula (1-10) which, of eom-se, becomes a better 
approximation the larger L becomes. If X = 5, ^ ^ 10, and S = 0.02, 
the smallest value, a, such that 

(X - a){n ~ a) = Xixe-'^'^'-"'' 

is a = 1.548. The asymptotic formula for F(L) now becomes 

F(L) ^ l.Oiae"''''^ 

which is found in Fig. 8 as ?/2 . 

S.L3 Bounds Using the Asymptotic Exponent 

Formulas (1-12) through (1-18) give a scheme for computing both 
upper and lower bounds for F{L) which have the right behavior for 
large L, and also agree with the solution at L = 0. They become 

FiD ^ l.OGTe"''*''' - 0.007e"'''', 

and 

F(L) ^ l.lOSc"''*'^ - O.lOSe-"'', 
respectively, and are represented by 1/3 and 1/4 in Figiu-e 8. 

5.1.4 ^n Upper Bound by a Discrete Markov Process 

If we mark on the positive .r-axis the points n8/2, n = 0, 1, 2, ■ ■ ■ , 
Ave can assign to each interval of length 5/2 thus created a state {ij), i, 
j = OY 1, as follows: i = if no point of the X-process is pre-sent in the 
interval, i = 1 if one or more points of the X-process are present, and 
similarly for j and m- An interval of length 8, made up of two adjacent 
intervals of length 5/2, may then be represented by a number between 
and 15 in binary notation, where 3, 6, 7, 9, and 11-15 represent a 
coincidence within the interval of length 5. We now define a Markov 
process as follows: in the interval S t < 8, let p/"', i = 0, 1, 2, 4, 5, 
8, 10, be the probabilities of occurrence of the z"" state, so that, for exam- 
ple, pr = .-"V^"', and pr = e-^V^d - e""'). These are the 
states in which there is no coincidence in (0, 8). In addition, let g"" 
represent the probability of all the other states put together; i.e., of a 
comcidence in (0, 5). We now define p/"', i = 0, 1, 2, 4, 5, 8, 10 as the 
probability of the z"" state in the interval {nS/2, (n + 2)5/2), where we 
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require in addition that all states in the intervals {l-h/2, (/,' + 2)5/2), 
/,■ < n, are from the same "no coincidence" index set. We define q " as 
the probabiUty of a state :^, G, 7, 0, or 11-15, in some interval (A-5/2, 
(k + 2)5/2), k ^ ?i.. There are then transition pioliabilities from states 
in the n - i'* to states in the n^ interval. For example, 



W"' = « 



KPo 



= e e 






), 



and 



^(.) ^ ^(.-1) ^ ^1 _ -*)(i _ ,-.«)(p, 



+ P4^"-^' + 



(n-l) 



+ (1 - e-")(p/"-» + Pa'"-") + {1 - C-''')(P2*"^" + Pio^""^'). 

The quantity 1 - g"" is then an upper bound for the probabiUty of no 
coincidences (upper because it is possible for a coincidence to occur in 
the process which is not counted in this subdi^'ision of it). The curve 
yz in Fig. 8 is drawTi through points at L =- ?i5/2 computed in this 
manner. 

To summarize the results, we see that the asymptotic formula and 
the lower bound are both indistinguishable from the right answer; the 
upper bounds are fairly far off. The upper bound derived by the Marko\- 
process is better than that derived from the integral equation until 
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Fig. 9 — Probability of no coincidences in a 26 X ^S square; neighborhoods are 
square. 
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about L = 0.5 (25 iterations), when the integral equation upper bound 
becomes better. 

5.2 A Single Pattern in a Square 

To tost our higher-dimensional bounds, we consider again coincidences 
in a single Poisson pattern in a square of side 25. The exact probability 
of no coincidences was given in Part IV assuming square neighborhoods. 
The lower bound (Sec. 4.2) 

applies using V = (26)' and S{5) = (25)" for square neighborhoods. To 
use the lower bound 1 - E we note that the exact expected number of 
coincidences is 

^ = 1'^' I I Mx,y)dxdy 

where A{x, y) is the area of the intersection of the given square with 
the square neighborhood centered at {x, y). The lower bound lal — E = 
1 - 9\V/2. The upper bound p*^ can be used if the square is cut into 
/C = 4 squares of side 5, each with a probabihty p - (1 -|- X5^) exp - 
X3 of no coincidence. 

These bounds, together with the exact probability, are plotted as 
functions of X5" in Fig. 9. When X5" is small, the 1 - _E bound is correct 
to terms of order 0(X''5^). This might have been predicted from (4-6) 
since it seems reasonable that Qi , Qi , ■■ ■ should be of higher order in 
X than Qi when X is small. Ultimately the first lower bound becomes a 
better estimate. It must be recognized that this other lower bound is 
being tested mider very severe conditions. Since every point of the 
square has a neighborhood which intersects the boundary, the errors 
from source (b) of Part V are considerable. 

The authors ^ish to thank D. W. Hagelbarger and H. T. O'Neil for 
their assistance in the course of the calculations reported in this section, 
and Miss D. T. Angell for preparing .some of the figures. 
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